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(N 
> 

We study the production of electroweak gauge bosons at high energies in the framework 
of fc^-factorization QCD approach. The amplitude for production of a single W ± or Z° boson 
associated with quark pair in the fusion of two off-shell gluons is calculated. Contributions 
i/-) . from the valence quarks are calculated using the quark-gluon interaction and quark-antiquark 
annihilation QCD subprocesses. The total and differential cross sections (as a function of the 
transverse momentum and rapidity) are presented and the ratio of cross sections for and 
Z° boson production is investigated. The conservative error analysis is performed. In the 
numerical calculations two different sets of unintegrated gluon distributions in the proton 
are used: the one obtained from Ciafaloni-Catani-Fiorani-Marchesini evolution equation 
and the other from Kimber-Martin-Ryskin prescription. Theoretical results are compared 
with experimental data taken by the D0 and CDF collaborations at the Tevatron. We 
demonstrate the importance of the quark component in parton evolution in description of 
the experimental data. This component is very significant also at the LHC energies. 



PACS number(s): 12.38.Bx, 12.15.Ji 



1 Introduction 

The theoretical and experimental studying the vector {W^ and Z°) boson production 
at high energies provide information about the nature of both the underlying electroweak 
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interaction and the effects of Quantum Chromodynamics (QCD). In many respects these 
processes have become one of most important "standard candles" in experimental high en- 
ergy physics [1-9]. At the Tevatron, measurements of W ± and Z° inclusive cross sections 
are routinely used to validate detector and trigger perfomance and stability. Data from 
gauge boson production also provide bounds on parametrizations used to describe the non- 
perturbative regime of QCD processes. At the LHC, such measurements can serve as a useful 
tool to determine the integrated luminosity and can also be used to normalize measurements 
of other production cross sections (for example, cross section of W + n-jets or diboson pro- 
duction). Additionally, studying of inclusive vector boson production is necessary starting 
point for investigations of Higgs or top quark production where many signatures can include 
these bosons. 

At leading order (LO) of QCD, W ± and Z° bosons are produced via quark- ant iquark 
annihilation. Beyond the LO Born process, vector boson can also be produced by q + g inter- 
actions, so both the quark and gluon distribution functions of the proton play an important 
role. Theoretical calculations of the W ± and Z° production cross sections have been carried 
out at next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) [10-14] of 
QCD. The NLO cross section is ~ 25% larger than the Born-level cross section, and the 
NNLO cross section is an additional ~ 3% higher. However, these perturbative calculations 
are reliable at high p T only since diverge in the small p T <C m region with terms proportional 
to Inm/pT (appearing due to soft and collinear gluon emission). Therefore, the soft gluon 
resummation technique [15-19] should be used to make QCD predictions at low p?. The 
traditional calculations combine fixed-order perturbation theory with analytic resummation 
and some matching criterion. The analytic resummation can be performed either in the 
transverse momentum space [20] or in the Fourier conjugate impact parameter space [21]. 
Differences between the two formalisms are discussed in [22]. 

An alternative description can be provided by the /^-factorization approach of QCD [23, 
24]. This approach is based on the familiar Balitsky-Fadin-Kuraev-Lipatov (BFKL) [25] or 
Catani-Ciafaloni-Fiorani-Marchesini (CCFM) [26] gluon evolution equations and takes into 
account the large logarithmic terms proportional to Inl/x. This contrasts with the usual 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [27] strategy where only the large log- 
arithmic terms proportional to ln/i 2 are taken into account. The basic dynamical quantity 
of the /^-factorization approach is the unintegrated (i.e., ky-dependent) parton distribution 
f a (x, ky, /i 2 ) which determines the probability to find a type a parton carrying the longitu- 
dinal momentum fraction x and the transverse momentum k T at the probing scale /z 2 . In 
this approach, since each incoming parton carries its own nonzero transverse momentum, 
the Born- level subprocess q + q' — > W ± /Z° already generate the pt distribution of produced 
vector boson. Similar to DGLAP, to calculate the cross sections of any physical process 
the unintegrated parton density /„( ) has to be convoluted [23, 24] with the relevant 

partonic cross section which has to be taken off mass shell (k^-dependent). The soft gluon 
resummation formulas are the result of the approximate treatment of the solutions of the 
CCFM evolution equation [28]. Other important properties of the /c^-factorization formal- 
ism are the additional contribution to the cross sections due to the integration over the k 2 - 
region above yU 2 and the broadening of the transverse momentum distributions due to extra 
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transverse momentum of the colliding partonsj. 

The fey-factorization formalism has been already applied [30] to calculate transverse mo- 
mentum distribution of the inclusive and Z° production at Tevatron. The calcula- 
tions [30] were based on the usual (on-mass shell) matrix elements of the quark- ant iquark 
annihilation subprocess q + q' — > W ± /Z° which embedded in precise off-shell kinematics. 
However, an important component of the calculations [30] is the unintegrated quark distri- 
bution in a proton. At present these distributions are only available in the framework of the 
Kimber-Martin-Ryskin (KMR) approach [31] since there are some theoretical difficulties in 
obtaining the quark densities immediately from CCFM or BFKL equation^ (see, for exam- 
ple, review [29] for more details). As a result the dependence of the calculated cross sections 
on the non-collinear evolution scheme has not been investigated. This dependence in gen- 
eral can be significant and it is a special subject of study in the fer-factorization formalism. 
Therefore in the present paper we will try a different and more systematic way. Instead of us- 
ing the unintegrated quark distributions and the corresponding quark-antiquark annihilation 
cross section we calculate off-shell matrix element of the g* + g* —>■ W ± /Z° + qq' subprocess 
and then operate in terms of the unintegrated gluon densities only. In this scenario, at the 
price of considering the 2^3 rather than 2 — > 1 matrix elements, the problem of unknown 
unintegrated quark distributions will reduced to the problem of gluon distributions. How- 
ever, since the gluons are only responsible for the appearance of the sea but not valence 
quarks, the contribution from the valence quarks should be calculated separately. Having in 
mind that the valence quarks are only important at large x, where the traditional DGLAP 
evolution is accurate and reliable, this contribution can be taken into account within the 
usual collinear scheme based on the q + g* — > W^/Z + q' and q + q' — > W^/Z matrix ele- 
ments convoluted with the on-shell valence quark and/or off-shell gluon densities^). Thus, the 
proposed way enables us with making comparisons between the different parton evolution 
schemes and parametrizations of parton densities^. 

We should mention, of course, that this idea can only work well if the sea quarks appear 
from the last step of the gluon evolution — then we can absorb this last step of the gluon 
ladder into hard matrix element. However, this method does not apply to the quarks coming 
from the earlier steps of the evolution (i.e., from the second-to-last, third-to-last and other 
gluon splittings). But it is not evident in advance, whether the last gluon splitting dominates 
or not. The goal of our study is to clarify this point. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas of 
the /^-factorization approach with a brief review of calculation steps and the unintegrated 
parton densities used. We will concentrate mainly on the g* + g* — > W ± /Z° + qq' subprocess. 
The evaluation of q + g* — > W ± /Z° + q' and q + q' — > W ± /Z° contributions is straightforward 
and, for the reader's convenience, we only collect the main relevant formulas in Appendix. In 
Section 3 we present the numerical results of our calculations. The central point is discussing 
the role of each contribution mentioned above to the cross section. Special attention is put 
on the transverse momentum distributions of the and Z° boson measured by the D0 [5, 
8, 9] and CDF [4] collaborations. Section 4 contains our conclusions. 

^Tor an introduction to /^-factorization, see, for example, review [29]. 
2 Unintegrated quark density was considered recently in [32]. 

3 To avoid the double counting we have not considered here q + <f — > W / Z° + g subprocess. 

4 The similar scenario has been applied recently to the prompt photon hadroproduction at Tevatron [33]. 
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2 Theoretical framework 



As the off-shell gluon-gluon fusion g* + g* — > W ± /Z° + qq' is calculated for the first time 
in the literature, we find it reasonable to explain it in more detail. 



2.1 Kinematics 

We start from the kinematics (see Fig. 1). Let p^ and p^ be the four- momenta of the 
incoming protons and p the four-momentum of the produced W ± /Z° boson. The initial off- 
shell gluons have the four-momenta k\ and k 2 and the final quark q and antiquark q' have the 
four- momenta p\ and p 2 and the masses m± and m 2 , respectively. In the pp center-of-mass 
frame we can write 

P {1) = ^(1,0,0,1), ^ = ^(1,0,0,-1), (1) 

where y/s is the total energy of the process under consideration and we neglect the masses 
of the incoming protons. The initial gluon four-momenta in the high energy limit can be 
written as 

h = x lP {1) + k 1T , k 2 = x 2 p {2) + k 2T , (2) 

where kyr and k 2 x are the transverse four- momenta. It is important that kf T = —k\ T ^ 
and k| T = — k 2T ^ 0. From the conservation laws we can obtain the following relations: 

kiT + k 2 T = Pit + P2T + Pt, 

xi\/s = miTe yi + m 2 Te y2 + m^e^, (3) 

x 2 \/s = m 1T e~ yi + m 2T e~ V2 + mxe~ y , 

where p?, itlt and y are the transverse momentum, transverse mass and center-of-mass 
rapidity of produced W ± / Z° boson, pyr and p 2 T are the transverse momenta of final quark 
q and antiquark q', y 1 , y 2 , m XT and m 2T are their rapidities and transverse masses, i.e. 
m\ T = mf + p? T . 



2.2 Off-shell amplitude of the g* + g* -> W ± /Z° + qq' subprocess 

There are eight Feynman diagrams (see Fig. 2) which describe the partonic subprocess 
g* + g* — > W ± /Z° + qq' at aa 2 s order. Let ei, e 2 and e be the initial gluon and produced 
gauge boson polarization vectors, respectively, and a and b the eight-fold color indices of the 
off-shell gluons. Then the relevant matrix element can be presented as follows: 

M 1 = 9 2 u(Pi)t a l^ P r^ 1+r ?\ 2 T^ z e x * 2 ~ff + m ° tYe^fo), (4) 

mf - (pi - h) 2 ' mi - (k 2 - p 2 y 

KA 2-, w6 v Pi-k 2 + m 1 x h-p 2 + m 2 

M 3 = / «fe ) t°7% " *' + ?' tS% 77^ + ""A £A „fe) , (6) 

mf — (p 1 — ki) z mf — {—p 2 — p) z 
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ka 2-/ \ + b v Pi-k 2 + m 1 -p 2 -p + m 1 A 

mj-(p 1 -k 2 ) 2 mj-(-p 2 -p) 2 

AW ufo) r^e, ^ _ ^ +p)2 * 7 e, mi _ {h _ p y 7 %.(p 2 ), (8) 

2 — / \ ttA Pl+P + m 2 , a „ k 2 -p2+m 2 , b v , s , nN 



-Mr = <? 2 «(pi) ^"(fei, fc 2 , -*i - h) - , 2 r fcc t c x 



(fci + fc : 



2y 



m( — (—p 2 — p) 1 



(10) 



K A 2-1 \rn\ yl ~ V ~ 1 10 2 

M 8 = g «(Pi)T w , z ex rn 2_ {pi+p)2 x (n) 

xyC^(h, k 2 , -h - k 2 )-^—f ab H c u{p 2 ). 

(h + k 2 y 

In the above expressions C^ p (k,p, q) and T^ z are related to the standard QCD three-gluon 
coupling and the W^ ± /Z°-fermion vertexes: 

C^(k,p, q) = gTip ~ k) p + 9 up (q ~ p f + <T (* ~ ?)", (12) 

T w= m e . « 7 A (l-7 5 )^, (13) 
2V2 smO w 



T z = —^T^ - 7 5 ) - 2e q sin 2 6 W \ , (14) 

sin zt>w L J 



where 1^ and e q are the weak isospin and the fractional electric charge (in the positron charge 
e units) of final-state quark q, 9 W is the Weinberg mixing angle and V qQ > is the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element. Of course, in the case of Z° production m\ 
equals m 2 . The summation on the W ± /Z° polarization is carried out by the covariant 
formula 

£e^(p)e*>) = -<r + ^. (15) 
In the case of the initial off-shell gluon we use the BFKL prescription [23, 24]: 

£^ K ^ = W (16) 

This formula converges to the usual expression v — —g pv after azimuthal angle aver- 

aging in the kr — * limit. The evaluation of the traces in (4) — (11) was done using the 
algebraic manipulation system Form [34]. We would like to mention here that the usual 
method of squaring of (4) — (11) results in enormously long output. This technical problem 
was solved by applying the method of orthogonal amplitudes [35]. 

The gauge invariance of the matrix element is a subject of special attention in the kr- 
factorization approach. Strictly speaking, the diagrams shown in Fig. 2 are insufficient 
and have to be accompanied with the graphs involving direct gluon exchange between the 
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protons (these protons are not shown in Fig. 2). These graphs are necessary to maintain the 
gauge invariance. However, they violate the factorization since they cannot be represented 
as a convolution of the gluon-gluon fusion matrix element with unintegrated gluon density 
The solution pointed out in [24] refers to the fact that, within the particular gauge (16), 
the contribution from these unfactorizable diagrams vanish, and one has to only take into 
account the graphs depicted in Fig. 2. We have successfully tested the gauge invariance of 
the matrix element (4) — (11) numericallj@. 



2.3 Cross section for the inclusive W ± /Z° production 

According to the fc^-factorization theorem, the inclusive W^/ Z° production cross section 
via two off-shell gluon fusion can be written as a convolution 



a( P + W±/Z° +*) = £/ ^V s (*i, k 2 1T , » 2 )dk 2 lT ^x 



x 



(17) 

/ ^f g (x 2 y 2T ,^)dklAa(g* + g*^ W±/Z° + qq>), 

J Xo ' ZTX 



where a(g*+g* — > W ± /Z°+qq') is the partonic cross section, f g (x, k^, /x 2 ) is the unintegrated 
gluon distribution in a proton and 0i and (p 2 are the azimuthal angles of the incoming gluons. 
The multiparticle phase space Ud 3 pi/2EiS^{J2p m ~ Y,P° nt ) is parametrized in terms of 
transverse momenta, rapidities and azimuthal angles: 



d 3 Pi 



7T 



2E, 2 rix a 2vr 
Using the expressions (17) and (18) we obtain the master formula: 



dPiT d Vi-^-- ( 18 ) 



a{p + p ^ W ±/ Z o + X) = Y;j * s)2 \M(g* + g*^ W±/Z + gg')l 2 x 
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x/ 5 (xi, k 2 T , /J, 2 )f g (x 2 , k 2 T , /j, 2 )dk 2 1T dk 2 2T dp 2 1T p 2 T dydy 1 dy 2 



2 ,.2w - ,. 1,-' ..2,,n/2 ,,_2 ,,2 ,,,„,„ ,,„ #t #2 #1 #2 

2tt 2tt 2tt 2tt ' 



(19) 



where |-M(g* + g* — > W^/Z + gg')| 2 is the off-mass shell matrix element squared and 
averaged over initial gluon polarizations and colors, ip\ and ip 2 are the azimuthal angles of 
the final state quark and antiquark, respectively. We would like to point out again that 
\M(g* + g* — > W^/Z + qq')\ 2 strongly depends on the nonzero transverse momenta k 2 r 
and k 2T . If we average the expression (19) over 0i and <p 2 and take the limit k 2 r — > and 
k 2T — > 0, then we recover the expression for the W ± /Z° production cross section in the 
collinear aa 2 approximation. 

The multidimensional integration in (19) has been performed by means of the Monte 
Carlo technique, using the routine Vegas [36]. The full C++ code is available from the 
authors upon request^. 



5 At the preliminary stage of the work we have made a cross-check of the matrix elements which have 
been calculated independently by M. Deak and F. Schwennsen. 
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lipatov@theory.sinp.msu.ru 
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2.4 The KMR unintegrated parton distributions 

In the present paper we have tried two different sets of unintegrated parton densities in 
a proton. First of them is the Kimber-Martin-Ryskin set. 

The KMR approach [31] is the formalism to construct parton distributions f a ( 
unintegrated over the parton transverse momenta k 2 -, from the known conventional parton 
distributions xa(x, /i 2 ), where a = g or a = q. This formalism is valid for a proton as well as 
a photon and can embody both DGLAP and BFKL contributions. It also accounts for the 
angular ordering which comes from coherence effects in gluon emission. The key observation 
here is that the \i dependence of the unintegrated parton distributions f a ( ) enters at 

the last step of the evolution, and therefore single scale evolution equations (pure DGLAP) 
can be used up to this step. In this approximation, the unintegrated quark and gluon 
distributions are given by the expressions [31] 




where P a b(z) are the usual unregulated leading order DGLAP splitting functions, and q(x, /i 2 ) 
and g(x,n 2 ) are the conventional quark and gluon densitieCl- The theta functions which 
appear in (20) and (21) imply the angular-ordering constraint A = jjl/ (yU+|k T |) specifically to 
the last evolution step to regulate the soft gluon singularities. For other evolution steps, the 
strong ordering in transverse momentum within the DGLAP equations automatically ensures 
angular ordering. It is important that the parton distributions f a ( ) extended now 

into the k 2 -, > jj? region. This fact is in the clear contrast with the usual DGLAP evolutional. 

The virtual (loop) contributions may be resummed to all orders by the quark and gluon 
Sudakov form factors 




Numerically, we have used the standard GRV (LO) parametrizations [37]. 

8 We would like to note that cut-off A can be taken A = |ky|//i also [31]. In this case the unintegrated 
parton distributions given by (20) — (21) vanish for k^ > fj, 2 in accordance with the DGLAP strong ordering 
in kj,. 



7 



where z max = 1 — z min = ///(//+ |pt|)- The form factors T a (ky,/i 2 ) give the probability 
of evolving from a scale k\ to a scale fi 2 without parton emission. In according with (22) 
and (23) T a (kf.,/i 2 ) = 1 in the k| > /j 2 region. 

Note that such definition of the f a (x, ky, /z 2 ) is correct for kf~ > /Xq only, where fi ~ 1 
GeV is the minimum scale for which DGLAP evolution of the collinear parton densities 
is valid. Everywhere in our numerical calculations we set the starting scale yU to be equal 
/i — 1 GeV. Since the starting point of this derivation is the leading order DGLAP equations, 
the unintegrated parton distributions must satisfy the normalisation condition 

a(x,fi 2 ) = J f a (x,k 2 T ,fi 2 )dk 2 T . (24) 
o 

This relation will be exactly satisfied if one define [31] 

f a (x,k 2 T^ 2 )\^ T <^ =a(x,/i 2 )T a (/i 2 ,/i 2 ). (25) 

2.5 The CCFM unintegrated gluon distribution 

The CCFM gluon density has been obtained [38] from the numerical solution of the 
CCFM equation. The function f g ( ) is determined by a convolution of the non- 

perturbative starting distribution fjj°\x) and the CCFM evolution kernel denoted by A(x, k 2 ., /i 2 ): 

f B (x, k 2 , /i 2 ) = / d 4ff\x')A(-^A . (26) 

J OC \X / 

In the perturbative evolution the gluon splitting function P gg (z) including nonsingular terms 
is applied, as it was described in [39] . The input parameters in (x) were fitted to reproduce 
the proton structure functions F 2 (x,Q 2 ). An acceptable fit to the measured F 2 values was 
obtained [38] with x 2 /ndf = 1.83 using statistical and uncorrelated systematic uncertainties 
(compare to ^ jndf ~ 1.5 in the collinear approach at NLO). 

3 Numerical results 

We are now in a position to present our numerical results. First we describe the theoretical 
uncertainties of our consideration. 

Except the unintegrated parton distributions in a proton f q ( ), there are several 

parameters which determined the overall normalization factor of the calculated W ± /Z° cross 
sections: the quark masses m\ and m 2 and the factorization and renormalization scales fip 
and jiR (the first of them is related to the evolution of the parton distributions, the other 
is responsible for the strong coupling constant). In the numerical calculations the masses 
of light quarks were set to be equal to m u = 4.5 MeV, rrid = 8.5 MeV, m s = 155 MeV 
and the charmed quark mass was set to m c = 1.5 GeV. We have checked that uncertainties 
which come from these quantities are negligible in comparison to the uncertainties connected 
with the scale and/or the unintegrated parton densities. As it is often done, we choose the 
renormalization and factorization scales to be equal: jiR = hf = A* = £ m T (transverse 
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mass of the produced vector boson). In order to investigate the scale dependence of our 
results we vary the scale parameter £ between 1/2 and 2 about the default value £ = 1. 
For completeness, we set mw = 80.403 GeV, mz = 91.1876 GeV, sin 2 9 W = 0.23122 and 
use the LO formula for the strong coupling constant a s (fi 2 ) with rif = 4 active quark 
flavors at Aqcd = 200 MeV (so that ot a (M|) = 0.1232). Note that we use a special choice 
Aqcd = 130 MeV in the case of CCFM gluon (a,(M§) = 0.1187), as it was originally 
proposed in [38]. 

Before we proceed to the numerical results, we would like to comment on the effect of the 
higher order QCD contributions [30] . It is well-known that the leading-order fc^-factorization 
approach naturally includes a large part of therr|§. It is a corrections which are kinematic 
in nature arising from the real parton emission during the evolution cascade. Another part 
of high-order contributions comes from the logarithmic loop corrections which have already 
been included in the Sudakov form factors (22) and (23). However, there are also the non- 
logarithmic loop corrections, arising, for example, from the gluon vertex corrections to Fig. 2. 
To take into account these contributions we will use the approach proposed in [30]. It was 
demonstrated that main part of the non-logarithmic loop corrections can be absorbed in the 
so-called K-factor given by the expression 



K{q + 4 -> W ± /Z°) ~ exp 



a s (fx 2 ) 2 
Ztt 



(27) 



where color factor Cf = 4/3. A particular choice fi 2 = p^/ 3 m 2//3 has been proposed [22, 
30] to eliminate sub-leading logarithmic terms. We choose this scale to evaluate the strong 
coupling constant a s (fi 2 ) in (27). 

We begin the discussion by presenting a comparison between the different contributions 
to the W^/Z total cross section. The solid, dashed and dotted histograms in Figs. 3 — 6 
represent the g*+g* — > W ± /Z° + qq', q+q* — > W ± /Z° + q' and q+q' — > W ± /Z° contributions 
to the rapidity distributions of gauge boson calculated at the Tevatron (Figs. 3 and 4) and 
LHC conditions (Figs. 5 and 6). It is important that in the last two subprocesses we take into 
account only the valence quarks within the usual collinear approximation. For illustration, 
we used here the KMR unintegrated gluon density. We found that the role of the gluon- 
gluon fusion subprocess is greatly increased at the LHC energy: it contributes only about 
2 or 3% of the valence quark component at the Tevatron and more than 40% at the LHC. 
Moreover, in the last case it dominates over the valence contributions at the central rapidities. 
The contribution of the valence quark-antiquark annihilation subprocess is important at the 
Tevatron and gives only a few percents at the LHC energy. As expected, the contribution of 
the q + g* — > W ± / Z° + q' subprocess is significant in the forward rapidity region, \y\ > 2. At 
this point, we can conclude that the gluon-gluon fusion becomes an important production 
mechanism at high energies and therefore should be taken into account in the calculations. 
However, we would like to note that there is an additional contribution which is not included 
in the simple decomposition scheme proposed above. As it was mentioned above, in this 
scheme it was assumed that sea quarks appear only at last gluon splitting and there is no 
contribution from the quarks coming from the earlier steps of the evolution (and we absorb 
last step of the gluon ladder into hard matrix element g* + g* — > W ± /Z° + qq'). It is 



3 See, for example, review [29] for more details. 
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not clear in advance, whether the last gluon splitting dominates or not. In order to model 
this additional component, we have repeated the calculations using the KMR unintegrated 
quark densities (20) and the quark- ant iquark annihilation q + <f — > W ± /Z° matrix element. 
But in these evaluations we omited the last term and keep only sea quark in first term 
of (20). Thus, we switch off the pure gluon component of the sea quark distributions and 
remove the valence quarks from the evolution ladder. In this way only the contributions 
to the f q (x, k. T ,fi 2 ) originating from the earlier (involving quarks) evolution steps are taken 
into account. So, the dash-dotted histograms in Figs. 3 — 6 represent the results of our 
calculations. We have found the significant (by about of 50%) enhancement of the cross 
sections at both the Tevatron and LHC conditions. Therefore in all calculations below we 
will consider this mechanism as an additional production one. Finally, taking into account 
all described above components, we can conclude that the gluon-gluon fusion contributes 
about ~ 1% to the total cross section at Tevatron and up to ~ 25% at LHC energies. 

Now we turn to the transverse momentum distributions of the W ± and Z° bosons. The 
experimental data for the transverse momentum distributions come from both the D0 [8, 
9] and CDF [4] collaborations at Tevatron. These data were obtained at center-of-mass 
energy <Js = 1800 GeV. Measurements were made for W — » Iv and Z — > l + l~ decays; 
so that we should multiply our theoretical predictions by the relevant branching fractions 
f(W — > Iv) and f(Z — > These branching fractions were set to f(W — - > Iv) = 0.1075 

and f(Z — > l + l~) = 0.03366 [40]. In Figs. 7 — 9 we display a comparison of the calculated 
differential cross sections da/dpr of the and Z° boson production with the experimental 
data [4, 8, 9] in the low pt region, namely pt < 20 GeV. Next, in Figs. 10 — 12, we demon- 
strate the and Z° transverse momentum distributions in the intermediate and high pt 
regions. Additionally, in Figs. 13 and 14, we plot the normalized differential cross section 
( 1 /a) da/dpT of the boson production. The solid and dashed histograms correspond 
to the results obtained with the CCFM and KMR unintegrated gluon densities, respec- 
tively. All contributions discussed above are taken into account. The dotted histograms 
were obtained using the quark-antiquark annihilation matrix element convoluted with the 
KMR unintegrated quark distributions in a proton (in this case the transverse mometum of 
the produced vector boson is defined by the transverse momenta of the incoming quarks). 
We found an increase in the cross section calculated in the proposed decomposition scheme 
(where only the unintegrated gluon densities used). In this scheme, we obtain that both 
the CCFM and KMR gluon distributions reproduce well the Tevatron data within the un- 
certainties, although the KMR gluon tends to slightly underestimate the data in the low 
Pt region. The difference between the solid and dashed histograms in Figs. 7 — 14 is due 
to different behaviour of the CCFM and KMR gluon densities. The predictions based of 
the quark-antiquark annihilation subprocess lie below the experimental data but agree with 
them in shape. This observation coincides with the one from [30] where an additional factor 
of about 1.2 was introduced to eliminate the visible disagreement between the data and 
theorj@. 

An additional possibility to distinguish the two calculation schemes comes from the study- 
ing of the ratio of the and Z° boson cross sections. In fact, since W and Z° production 

10 In Ref. [30] authors have explained the origin of this extra factor by the fact that the input parton 
densities (used to determine the unintegrated ones) should themselves be determined from data using the 
appropriate non-collinear formalism. 
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properties are very similar, as the transverse momentum of the vector boson becomes smaller, 
the radiative corrections affecting the individual distributions and the cross sections of hard 
process are factorized and canceled in this ratio. Therefore the results of calculation of this 
ratio in the decomposition scheme (where the 0(aa s ) and 0(aa 2 s ) subprocesses are taken 
into account) and the predictions based on the 0(a) quark- ant iquark annihilation should 
differ from each other at moderate and high values. This fact is clearly illustrated in 
Fig. 15 where the ratio of W ± and Z° cross sections as a function of the transverse momen- 
tum is displayed. As it was expected, there is practically no difference between all plotted 
histograms in the low p T region. 

As a final point of our study, we discuss the scale dependence of our results. In Figs. 16 
and 17 we show the total cross section of the W ± and Z° boson production as a function of 
the total center-of-mass energy y/s. Here, the solid and dotted histograms correspond to the 
results obtained with the CCFM and KMR unintegrated gluon densities, respectively. The 
the upper and lower dashed histograms correspond to the scale variations in the CCFM gluon 
density as it was described above. We find that the scale uncertainties are the same order 
approximately as the uncertainties coming from the unintegrated gluon distributions. This 
fact is the typical one for the leading-order /cy-factorization calculations. Our predictions 
for the W ± and Z Q boson total cross section agree well with the data in a wide range. 

4 Conclusions 

We have studied the production of electroweak gauge bosons in hadronic collisions at high 
energies in the /^-factorization approach of QCD. Our consideration is based on the scheme 
which provides solid theoretical grounds for adequately taking into account the effects of 
initial parton momentum. The central part of our derivation is the off-shell gluon-gluon 
fusion subprocess g* + g* — > W ± /Z° + qq'. At the price of considering the corresponding 
matrix element rather than q + q' — » W ± /Z° one, we have reduced the problem of unknown 
unintegrated quark distributions to the problem of gluon distributions. This way enables us 
with making comparisons between the different parton evolution schemes and parametriza- 
tions of parton densities. Since the gluons are only responsible for the appearance of sea, 
but not valence quarks, the contribution from the valence quarks has been calculated sep- 
arately. Having in mind that the valence quarks are only important at large x, where the 
traditional DGLAP evolution is accurate and reliable, we have calculated this contribution 
within the usual collinear scheme based on q + g* — > W ± /Z° + q' and q + q' — > W ± /Z° 
partonic subprocesses and on-shell parton densities. 

We have studied in detail the different production mechanisms of W ± and Z Q bosons. 
We find that the off-shell gluon-gluon fusion gives ~ 1% and ~ 25% contributions to the 
inclusive W^/Z production cross sections at the Tevatron and LHC. Specially we simulate 
the contribution from the quarks involved into the earlier steps of the evolution cascade (i.e., 
into the second-to- last, third-to-last and other gluon splittings) and find that these quarks 
play an important role at both the Tevatron and LHC energies. It was demonstrated that 
corresponding corrections should be taken into accout in the numerical calculations within 
the /c^-factorization approach. 

We have calculated the total and differential W ± and Z° production cross sections and 
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have made comparisons with the Tevatron data. In the numerical analysis we have used 
the unintegrated gluon densities obtained from the CCFM evolution equation and from the 
KMR prescription. Our numerical results agree well with the experimental data. 

When the present paper was ready for publication, we have learned about the results 
obtained by M. Deak and F. Schwennsen [41], who used the same theoretical approach, but 
focused attention on slightly different aspects of the problem. These authors concentrate on 
the associated W ± /Z° production with heavy quark pairs (mainly on the Zbb final state), 
where the gluon-gluon fusion subprocess dominates. In additional to that, we consider 
quark subprocesses, which are important for inclusive W ± /Z° production. We show that 
the experimental data can be described with taking quark contributions into account. 
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6 Appendix A 

Here we present the compact analytic expressions for the cross section of the W ± /Z° 
production via q + q' — > W ± /Z° subprocess in the /c^-factorization approach. Let us define 
the transverse momenta and azimuthal angles of incoming quark q and antiquark q' as 
k 1T and k 2T and 4>i and <p 2 , respectively. The produced vector boson has the transverse 
momentum p T (p T = k 1T + k 2T ) and center-of-mass rapidity y. The W ± /Z° production 
cross section can be written as 

a(p + p^W ± /Z° + X)=J2 I 7 ^ , 2 \M(q + q' ^W ± /Z°)\ 2 x 
xf q (x 1: k 2 1T , fi 2 )f q (x 2 , k 2 2T , fi 2 )dk 2 1T dk 2 2T dy d ^ , 

where f q (x, k^, /i 2 ) is the unintegrated quark distributions given by (20). In the high-energy 
limit the fractions X\ and x 2 of the initial proton's longitudinal momenta are given by 

Xiy/s = rriT e v , x 2 y/s = tjit e~ y . (^4-2) 

where m T is the transverse mass of the vector boson. The squared matrix element \Ai{q + 
q' — > W^)! 2 summed over final polarization states and averaged over initial ones is 



\M(q + q' - W±)\ 2 = - / , [(mi - m 2 ) 2 + m 2 w {m{ + m 2 ) - 2m i w \ , (A.3) 
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where m\ and m 2 are the masses of incoming quarks. In the case of Z° boson production, 
the squared matrix element \A4(q + q — > Z°)\ 2 summed over final polarization states and 
averaged over initial ones is 



\M(q + q^Z 



0\ 1 2 



2e 2 



9 sin 2 29 w 



■x 



x 



:m, 



m 2 ) [li q) ] + 2e q (2m 2 + m|) sin 2 9 W (e q sin 2 9 W - 



(A4) 



where m, e q and 



Jg 9 ^ is the mass, fractional electric charge and weak isospin of incoming 



quark. Note that there is no obvious dependence on the transverse momenta of the initial 
quark and antiquark. However, this dependence is present because the true off-shell kine- 
matics is used. In particular, the incident parton momentum fractions X\ and x 2 have some 
dependence. If we take the limit k 2 T — > and k 2 T — > 0, then we recover the relevant 
expression in the standard collinear approximation of QCD. 



7 Appendix B 

Here we present the analytic expressions for the cross section of the W ± /Z° production 
via q + g* — > W ± /Z° + q' subprocess in the /c T -factorization approach. Let us define the 
transverse momenta and azimuthal angles of incoming quark and off-shell gluon as and 
k 2 T and <\>\ and 02, respectively. In the following, s, i and u are usual Mandelstam variables 
for 2 — > 2 subprocess. The W ± /Z° production cross section can be written as follows: 
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x/ 9 (xi, k 2 T , jj 2 )f g (x 2 , k|r» i2 2 )dk 2 1T dk 2 2T dpTdydy'—± 



2tt 2tt ' 



where y' is the rapidity of the final quark g'. The fractions X\ and ^2 of the initial proton's 
longitudinal momenta are given by 

xi^/s = rriT e y + m' T e y , X2\/s = tut e~ y + m' T e~ y . {B.2) 

where ttit and m' T are the transverse masses of the vector boson and final quark q' . If 
we take the limit k 2 T — > and k| T — > 0, then we recover the relevant expression in the 
usual collinear approximation. The squared matrix elements \M(q + g* — > + g')| 2 an d 
|.M(g + g* — > Z° + <?)| 2 summed over final polarization states and averaged over initial ones 
are 

w , + ^^^ +rt |. = -^__^_ , (B ,) 
+ 9 - _ + = {BA) 

where 

F w = -8(m^(3m2 - i) + m\{m\ + m 2 2 {2m 2 w - 5s - 7t) + t(s + 2t)) + 



13 



m\{m\ + m\{8m 2 w - 3(s + m 2 2 (-Qm 4 w + 3s 2 + 13s* + 5t 2 - 6m 2 w (s + *))- 
i(-2m 4 v - 2m^s + (s + *) 2 )) + m 2 (3m^ + m 6 2 (2m 2 w -7s- 5*) + 
m 4 2 (-Qm 4 w + 5s 2 + 13s* + 3* 2 - 6m^(s + *)) + m 2 2 (4m 6 w - s 3 - lis 2 * - list 2 - i 3 + 

6m 4 w (s + *) + Qm 2 w (s 2 + t 2 )) + *(- 4 "4 + 2 <« + s(s + i) 2 - 
2m 2 w {2s 2 - st + t 2 ))) + s(-m 8 + m 6 (2s + *) + 2m 2 w i{2m 4 v + s 2 + t 2 - 2m^(s + *))+ 
m\(2m 4 w + 2m^* - (s + t) 2 ) + m 2 (-4m^ + 2m^* + *(s + i) 2 - 
2m 2 w (s 2 - st + 2t 2 ))) + (ml + m 8 2 + m\{m 2 w - 2(s + *)) + m\(m 2 w - 2(s + *))+ 

m 4 2 {-2m 4 w - 2m^* + (s + t) 2 ) + m 2 2 m 2 w {hs 2 + 4s* + P + m 2 w (-8s + 4*))- 
2m^(2s*(s + *) + m^(s 2 - 4si + t 2 )) + m\(-2m\ - 2m 4 v - 2m 2 w s + (s + *) 2 + 
m 2 (m^ + 2(s + *))) + m 2 (m 2 y (s 2 + 4m^(s - 2*) + 4st + 5t 2 ) + m^m^ + 2(s + *))+ 
m 2 2 (8m 4 w - 6m 2 w (s + *) - 2(s + t) 2 )))(-k 2 T ) + Am 2 w (m\ - s)(m 2 - *)k 4 T ), (B.5) 
F z = -2e q I ( 3 q) m 2 z sin 2 6 w {6m 8 - s*(2m| + s 2 + i 2 - 2m|(s + *))- 
m 4 (2m 4 z + 3s 2 + 14s* + 3t 2 - 2m|(s + + m 2 (s 3 - 8m|s* + 7s 2 i+ 
7 si 2 + i 3 + 2m|(s + *))) + 2e 2 m| sin 4 # w (6m 8 - s*(2m| + s 2 + t 2 - 2m 2 z (s + *))- 
m 4 (2m z + 3s 2 + 14s* + 3t 2 - 2m|(s + *)) + m 2 (s 3 - 8m 2 z si + 7s 2 i + 7s* 2 + i 3 + 



2m%(s + *))) + [4 9) J" (-4m 10 + m 8 (-6m| + 8(s + *)) - m 2 z si(2m% + s 2 + t 2 - 

2m 2 z (s + t)) + m 6 (Qm z - 5s 2 - Usi - 5t 2 + 6m|(s + £))+ 
m 4 (-2m| + s 3 + 7s 2 t + 7si 2 + £ 3 - 4m|(s + t) - m 2 z (3s 2 + 2st + 3t 2 )) + 

m 2 (-2m 4 z si + 2m%(s + t) - si(s + i) 2 + m|(s 3 + s 2 t + st 2 + f 3 ))) + 
m 2 z (2e q I^ q) sin 2 ^ w (2m 4 (m| - s - i) - 2si(s + t) - m 2 z (s 2 - Asi + i 2 )- 
2m 2 (-Asi + m 2 z (s + t))) + 2eJ sin 4 w (2s*(s + t) + 2m 4 (-m| + « + + 

m|(s 2 - 4s* + i 2 ) + 2m 2 (-Asi + m 2 z (s + £))) + [4 9) ] 2 (-2m 6 + 2si(s + t) + 
m 2 z {s 2 - 4s* + t 2 ) + m 4 (-2m| + 4(s + *)) + m 2 (-3s 2 - 4s* - 3* 2 + 
2m|(s + *))))(-k 2 T ) - 2m|(m 2 - s)([/^] 2 - 

2e g /i 9) sin 2 + 2e 2 sin 4 O w ) (m 2 - *)k 4 T . (B.Q) 
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Figure 1: Kinematics of the g* + g* — > W^/Z + qq' process. 
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Figure 2: Feynman diagrams which describe the partonic subprocess g* + g* — > W ± /Z° + qq' 
at the leading order in a s and a. 
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Figure 3: Different contributions to the inclusive W ± boson production at the Tevatron 
conditions. The solid, dashed and dotted histograms correspond to the g* + g* — > W ± + qq 1 , 
q + g* — > W ± + q' and q + q' — > W ± subprocesses. In the last two cases only the valence 
quarks are taken into account. The dash-dotted histogram represents contribution from the 
quarks coming from the earlier steps of the evolution. The thick solid histogram represents 
the sum of all contributions. The KMR unintegrated parton densities in a proton are used. 
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Figure 7: Transverse mometum distribution of the W ± boson production. The solid and 
dashed histograms correspond to the results obtained with the CCFM and KMR uninte- 
grated gluon densities in a proton, respectively. The dotted histograms were obtained by 
using the quark- ant iquark annihilation matrix element convoluted with the KMR uninte- 
grated quark distributions. The cross sections times branching fraction f(W — > lu) are 
shown. The experimental data are from D0 [9]. 
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Figure 8: Transverse mometum distribution of the Z Q boson production. Notations of 
histograms are the same as in Fig. 7. The cross sections times branching fraction f(Z — > 
l + l~) are shown. The experimental data are from D0 [8]. 
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Figure 9: Transverse mometum distribution of the Z Q boson production. Notations of 
histograms are the same as in Fig. 7. The cross sections times branching fraction f(Z — > 
l + l~) are shown. The experimental data are from CDF [4]. 
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Figure 10: Transverse mometum distribution of the W ± boson production. Notations of 
histograms are the same as in Fig. 7. The cross sections times branching fraction f(W — > lu) 
are shown. The experimental data are from D0 [9]. 
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Figure 11: Transverse mometum distribution of the Z Q boson production. Notations of 
histograms are the same as in Fig. 7. The cross sections times branching fraction f(Z — > 
l + l~) are shown. The experimental data are from D0 [8]. 
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Figure 12: Transverse mometum distribution of the Z Q boson production. Notations of 
histograms are the same as in Fig. 7. The cross sections times branching fraction f(Z — > 
l + l~) are shown. The experimental data are from CDF [4]. 
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Figure 14: Normalized transverse mometum distribution of the W ± boson production. No- 
tations of histograms are the same as in Fig. 7. The earlier experimental data are from 
D0 [5]. 
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Figure 15: Ratio of differential cross section for W ± to Z° production. Notations of his- 
tograms are the same as in Fig. 7. The experimental data are from D0 [8]. 
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Figure 16: The total cross section of the inclusive W ± boson production as a function of \fs. 
The solid and dotted histograms correspond to the results obtained with the CCFM and 
KMR unintegrated gluon densities, respectively. The upper and lower dashed histograms 
correspond to the scale variations in CCFM gluon density as it was described in text. The 
cross sections times branching fraction f(W — > Iv) are shown. The experimental data are 
from UA1 [1], UA2[2], D0 [6, 7] and CDF [3]. 
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Figure 17: The total cross section of the inclusive Z° boson production as a function of 
yt~s. Notations of histograms are the same as in Fig. 16. The cross sections times branching 
fraction f(Z — > l + l~) are shown. The experimental data are from UAl [1], UA2[2], D0 [6, 
7] and CDF [3]. 
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